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A rapidly rotating Bose-Einstein condensate in a symmetric two-dimensional trap can be described 
with the lowest Landau-level set of states. In this case, the condensate wave function xf}{x,y) is a 
Gaussian function of r^ = x'^ + y^, multiplied by an analytic function P{z) of the single complex 
variable z = x+iy; the zeros of P{z) denote the positions of the vortices. Here, a similar description is 
used for a rapidly rotating anisotropic two-dimensional trap with arbitrary anisotropy {uix/i^y < 1). 
The corresponding condensate wave function tp{x,y) has the form of a complex anisotropic Gaussian 
with a phase proportional to xy, multiplied by an analytic function P{C), where ^ (X x + i(5-y and 
< /3- < 1 is a real parameter that depends on the trap anisotropy and the rotation frequency. The 
zeros of P(C) again fix the locations of the vortices. Within the set of lowest Landau-level states 
at zero temperature, an anisotropic parabolic density profile provides an absolute minimum for the 
energy, with the vortex density decreasing slowly and anisotropically away from the trap center. 

PACS numbers: 03.75.Hh, 05.30.Jp, 67.40.Db 

I. INTRODUCTION 

The experimental creation of rapidly rotating Bose-Einstein condensates (BEC) generally involves anisotropic rotat- 
ing trap potentials P,S|3, yet most theoretical analyses of such systems have relied on an isotropic trap 0,|SIEI3)IS0- 
As emphasized by Ho the low-lying states in a symmetric two-dimensional trap are closely analogous to those in 
the lowest Landau level for a charged particle in a uniform magnetic field. This analogy allows a simplified description 
in the limit that the rotation frequency approaches the frequency loq of the symmetric confining trap. If the typical 
interaction energy is small compared to the spacing « 2hujQ between adjacent Landau levels, then the condensate 
wave function ■0 of the interacting rotating BEC can be constructed as a linear combination of the lowest Landau-level 
(LLL) states. The nth such state is simply proportional to z", where z = x + iy, multiplied by the ground-state 
Gaussian. It follows that such a linear combination involves an analytic function of z that is usually approximated by 
a polynomial P{z) oc X\j[z — Zj), where the zeros Zj of the polynomial represent the positions of the vortices in the 
two-dimensional condensate. If the vortex density is strictly uniform, then the overall density profile is also Gaussian 
with an effective condensate radius that grows and ultimately diverges as — > wq 0. In fact, this system can lower 
its energy by slightly reducing the vortex density near the outer edge of the condensate, and the actual density profile 
has a quadratic shape (an inverted parabola) [MIESISj as shown by both analytical and numerical studies. 

The c|ua ntum-mechanical problem of a particle in a rotating two-dimensional anisotropic harmonic trap is exactly 
soluble 10,'Tl', l^l, although the corresponding eigenstates have not been discussed previously in full detail. We here 
amplify Valatin's description |0 to construct the anisotropic analogs of the LLL states. Each such state Lpno{x,y) 
involves the anisotropic complex Gaussian ground-state eigenfunction (p()Q{x,y), multiplied by a polynomial Pn(C)j 
where C cx a; -I- iP-y is a single "stretched" complex variable, and < /3_ < 1 is a real parameter that depends on 
the trap anisotropy and the rotation frequency. Thus a linear combination of these LLL states for an anisotropic trap 
again involves an analytic function of the single complex variable C, (apart from the common overall factor ^pm). The 
corresponding zeros again represent the positions of the vortices, now in the rotating anisotropic BEC. 

Section II focuses on the eigenstates of the single-particle Hamiltonian Hq for a rotating anisotropic harmonic 
potential, starting with the classical trajectories and then obtaining the explicit form of the low-lying quantum 
mechanical states ipno that are the analogs of the lowest Landau-level states for a charged particle in a magnetic field. 
As in that case, the expectation value of both Hq and the angular momentum in the lowest Landau level can be 
reduced to corresponding expectation values of x"^ and y"^. The interacting dilute Bose-Einstein gas in this rapidly 
rotating anisotropic trap is treated in Sec. III. For a trial state constructed as a linear combination of ip„o{x,y), 
an anisotropic parabolic density profile provides the absolute minimum of the energy. The density of vortices is 
constant near the center but decreases slowly toward the edge of the condensate. Section IV contains a discussion 
and suggestions for additional study. 
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II. SINGLE-PARTICLE EIGENSTATES 



Consider a particle of mass m in an anisotropic two-dimensional harmonic potential (for definiteness, I assume 
oscillator frequencies tOx < ujy) that rotates uniformly at an angular velocity = perpendicular to the plane of 
the motion. In the rotating frame, this potential is time- independent, with the Hamiltonian 

Ho = + \m K + a;2 ^2) _ ^ ^^^^ _ ^^^^ ^ 

where the last factor involves the angular momentum — xpy — ypx- If < (^j,, the centrifugal force preferentially 
expands the condensate along the x axis. Although this case {lOx < ojy) is of principal interest here, it will also be 
valuable to see how the more familiar symmetric case emerges in the limit lOx = — ^o- 



A. Classical dynamical trajectories 

The normal modes of Eq. ^ are readily determined to have the frequencies 0, 0, 0, ITsL 



+n^^^\{ul-u:lY + ^u:ln\ (2) 

2 + '^y) mean-squared oscillator frequency. This general result contains several important 



1 (, ,2 

limits. 



1. In the symmetric case uox — — ^Oi the plus (minus) modes have frequencies uj± = T Specifically, 
the plus mode with frequency cj+ — uj^ ~ ^ has a reduced frequency when viewed from the rotating frame (as 
is evident physically) and a positive angular momentum (which explains the notation). Correspondingly, the 
minus mode has an increased frequency oj- — ujq + ^ and a negative angular momentum. 

2. li ujx < ^y, then the modes are nondegenerate even for $1 = 0, when they reduce to a;+ = lo^ and a;_ — ujy. 

3. For an anisotropic trap (uj^ < ^y) and rapid rotation with 5—1 — ^/ujx 0^, the plus normal- mode frequency 
vanishes, with 



. 2C.2 (c.2 

(jJ. 



2^ xyy x^^^ (3) 



3^2 -f ul 

Thus —duj+/dfl diverges for small S like 6"^^'^. In contrast, the minus normal-mode frequency remains finite at 
(5 = 0, with 

uj'^^SujI+ujI. (4) 

4. In the case of a rapidly rotating nearly symmetric trap with two small parameters |12| 6=1 — ^/uJx and 
77 = (jjy/ujx — Ij these eigenfrequencies simplify to 

"'+-^5^+^). —^2 + \^-5. (5) 



LOx OJ, 



X 



Note the sensitivity of in Eq. to the order of limiting procedures: (i) first 77 = (namely ujx = ujy = ujq) 
and then 5 — > or (ii) 5 ^ Q (namely — > ujx) at fixed Ux < For a symmetric trap, the plus frequency 
vanishes linearly with the small parameter 5] in contrast, for the anisotropic trap with uJx < the plus 
frequency vanishes like ^/S, with a coefficient proportional to y^. 

Figure 1(a) shows the two different normal-mode frequencies u}± (normalized to w^) as functions of ^/lOx for the 
typical anisotropy ujy/ojx = 1.2. 

For a symmetric trap with ojx = ^y ^ ^a, the plus (minus) normal modes with frequencies uj± = cjq T have 
counterclockwise (clockwise) circular orbits with positive (negative) helicity. To understand the physics of the two 
normal modes in the more general anisotropic case Ux < (^y, consider first the motion in the plus mode [lol |. It has 
the form 

x+{t) = xoe-'^+', y+{t) = 7/3+xoe-^"+*, (6) 
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FIG. 1: Behavior of relevant dimensionless quantities as a function of dimensionless rotation speed for the typical 

anisotropy ujy/uOx = 1-2. (a) dimensionless normal mode frequencies lo+/uJx and uj-/uJx and dimensionless auxiliary frequency 
"t/uLix\ (b) dimensionless polarization parameters /3+ and 13- . 



where < /3+ < 1 is a real non-negative dimensionless parameter and xq is an arbitrary amplitude. A detailed 
analysis shows that the parameter /3+ has two alternative representations 



In the rotating frame, the plus orbit is an ellipse with major axis oriented along x; the polarization parameter /3+ 
gives the ratio of the minor to major axes. The second form of Eq. {Tj) shows that uj+/(3+ « {ojy — uj'^)/{2uJx) for 

(5 ^ 0, so that u)+ and both vanish in this limit like VS. Physically, this behavior reflects the rotation-induced 
cancellation of the harmonic confinement in the x direction; the orbit then becomes linearly polarized as f2 — > Wa,. 
The plus motion is counterclockwise, with positive helicity and positive angular momentum. 
Similarly, the orbit for the minus mode has the parametric representation 

x^it) = z/3_yoe'^"-*, y-it) = yoe-'^-\ (8) 

where yo is an arbitrary amplitude and < /3_ < 1 is a real non-negative parameter with two alternative representa- 
tions ,1^ 

Unlike /3+, these relations show that has a nonzero limit /3_ w 2a;^/(3(jj^ + ij^lY^^ < 1 as (5 ^ 0. In the rotating 
frame, the minus orbit is an ellipse with major axis oriented along y; the polarization parameter /3_ gives the ratio 
of the minor to major axes. The minus motion is clockwise, with negative helicity and negative angular momentum. 
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It is again instructive to specialize these results to the case of a rapidly rotating nearly symmetric trap T?] , making 
use of the small parameters S — 1 — n/tOx and 77 — ajy/cux — 1. A straightforward expansion yields the approximate 
polarization parameters 



Although /3_ varies smoothly in this limit, the corresponding polarization for the positive mode has a singular 
behavior that depends on the relative magnitude of the two small parameters. More generally, /3+ and /?_ both have 
the value 1 in the limit of a symmetric trap, independent of Q, since the resulting motion is circularly polarized. 
Figure 1(b) shows the dependence of the two parameters /3± on fl/cux for an anisotropy ujy/cux = 1-2. Note the 
singular slope of uj+ and /3+ near the upper limit. 

B. Bogoliubov canonical transformation to diagonal Hamiltonian 

The structure of Hq in Eq. ^ is unusual because the term fiL^ = il {ypx — xpy) couples the otherwise independent 
X and y motions. This situation can be clarified by introducing the conventional ladder operators |15| 

1 / X dxPx \ t 1 f X dxPx \ 



where dx = \/'hl (muJx) , an d similarly for ay and a]^. With these operators, it is straightforward to see that the term 
riLz is proportional to [ij 

in [{uy + LOx) {alay - aj^a^,) + {ujy - lUx) {axUy - a^al)] . (12) 

The first term is "diagonal" in the creation and annihilation operators, but the second is "off-diagonal," similar to 
Bogoliubov's approximate Hamiltonian for a dilute Bose-Einstein gas [ig. Unfortunately, a direct diagonalization 
based on these "particle" operators involves considerable algebraic complexity pjL . 

Thus, it is preferable to return to the original single-particle Hamiltonian in Eq. fQ. Since Hq is quadratic in the 
coordinates and momenta, it can be diagonalized with a canonical transformation to new variables that obey the 
same Poisson brackets (in the classical case) or the same commutators (in the quantum case). Specifically, I follow 
Valatin |Tol | and introduce the generating function 

S{x, y- Q+, Q_) = -TO7 [A+A_Q+Q_ + \{^\ + A^) xy - X+Q+y - X-Q-x] . (13) 

Here, Q± are new canonical coordinates, X± are dimensionless constants given by 

^ cj± ^ .^^ith ^± = cj± + /3+/3_cjzp, (14) 

and 7 = ^+/f3+ — fi-//3- has the dimensions of a frequency. It follows from Eqs. (O, 0, © and H14(l that 7 has 
various equivalent representations 

7 = — +UJ-I3- ^ — +UJ+I3+ 
P+ P- 



20 n • (1^) 

For a symmetric trap with lOx — ^y = ^oj this frequency reduces to 7 = 2a;o for all O. In contrast, for an anisotropic 
trap, 7 diverges as O becomes small and approaches the value « (3a;^ -I- w^)/(2w2,) > 2^2, for n ujx- Figure 1(a) 
shows the normalized parameter ^ jwx as a function of n/ojx for ujy/ujx — 1.2. 

According to the general theory of classical Hamiltonian dynamics [T^ . any function like S{x,y;Q^,Q-) that 
depends on both the old and new coordinates will automatically generate a canonical transformation from old canonical 
variables to new canonical variables, with the corresponding momentum variables given by 

_ as _dS_ 

P - P - 9^ (16) 
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The first set of equations immediately yields the relations 

'Xl + Xl\ _ py 



2A+ / mjX+ ' 

that express the new coordinates Q± as linear combinations of the original coordinates and momenta. Similarly, the 
second set of equations can be used to find the corresponding relations for the new momenta P± 

P+ = TO7A+ ^ „ " - 1 I y + X+px, 



' Xl + X i 

For future reference, note the following alternative relations 



P = m-/X^ { l]x + X-Py. (18) 



P_ 

X = X+Q- 



they express the original coordinates in terms of the new canonical variables and will be valuable in the subsequent 
analysis. 

It is now straightforward to verify that the Hamiltonian has the following simple diagonal form when expressed in 
the new canonical variables 

H,^g- + lm.lQl + ^ + lrn.lQl. (20) 

One strategy is to substitute Eqs. H17|l and 1)18(1 directly into H20|l . which eventually reproduces the original Eq. Q. 
This new Hamiltonian 1(20(1 has the great advantage of immediately providing a quantum description of two inde- 
pendent harmonic oscillators with mass m and frequencies lo±. The corresponding quantum-mechanical annihilation 
operators a± and creation operators a\. follow from general quantum theory 

1 fQ±^.d±P±\ t 1 fQ± ■d±P±\ 



where d± = y^h/{muj±) are the oscillator lengths for the two separate modes. These operators obey the usual 

commutation relations [a±,Q!lj_] = 1 (all other commutators vanish). Note that d+ diverges as ^ cux, whereas 
remains finite in the same limit. In terms of these operators, the Hamiltonian takes the form |ll| 

Ho = \huj+ {a\a+ + a+a^^ + \fiuj- {a^_a^ + a^a^_^ . (22) 

C. Lowest Landau-level single-particle states for rotating anisotropic trap 

The single-particle ground state lpqq is given by the prescription a±LpQQ — 0, which leads to the explicit representation 



^oocxexp(-l^-l^) (23) 



as a Gaussian function of the two new coordinates Q± . Correspondingly, the complete set of normalized single-particle 
states ipn^n- is specified by two non-negative integers n± 
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Here, the eigenstate state ^Pn^n^ has {n^) quanta with frequency (w-); the energy eigenvalue is 

= ffjJ+ {n+ + 5) + Tiio- {n- + i) . (25) 
This eigenstate (pn+n^ has an angular momentum [Tll | 



Since duj^/dH. is negative {doj^/dft is positive), this result confirms that the plus (minus) mode has positive (negative) 
angular momentum. For the special case of a symmetric trap, the detailed form of these eigenstates is well known |l8l 

It is important to re-express the ground-state wave function ipoo in terms of the original canonical coordinates x and 
y. Valatin ^3 uses the generating function S{x,y;Q+,Q-) in Eq. (|13|l to obtain the explicit (factorized) expression 



ipoo{x,y) cx exp 



mj {P+x^ + P-y^ 



mxy 
exp I — - — 



7 If LUj^ LU 



2n(l + /3+/3_) 

Each of these two factors has an interesting structure. 

1. The first factor is a real anisotropic Gaussian with characteristic lengths ax and ay given by 



(27) 



2 _ 1 -I- h 2 _ 1 + P+P- ^ 



K - """" — , = ■ (28) 

In the limit of a symmetric trap, this Gaussian ground state becomes (poo{x,y) oc exp [— mwo {^^ + U^) /(2^)] , 
with the expected oscillator length do = y/Ti/ (towq)- More generally, the ground-state density for a rotating 
anisotropic trap is an anisotropic Gaussian, given by the corresponding normalized wave function 

1 f x^ y^ 



Voo{x, y) I = exp -] . (29) 



For a rapidly rotating anisotropic trap [lOx < ijJy and 5 ~ I — VL/lOx —f 0), the length ax diverges because /3+ 0, 
but ay remains finite. In this limit, the ground-state density becomes an essentially one-dimensional strip with 
Gaussian transverse profile and finite width ay 20, 21j. 

2. The second factor of lpqq involves a complex phase proportional to xy, which reflects the irrotational flow induced 
by the rotating anisotropic trap [2^ I23I l23|. The factor in square brackets (the coefficient of imxy/fi) has a 
rather intricate structure. It vanishes for a symmetric trap, because I3± = 1 and lo+ + iO- = 7 = Iujq. It also 
vanishes for a stationary anisotropic trap, but this limit requires a detailed analysis because each term separately 
diverges as $7 ^ 0. This phase will be seen to play an essential role in the following construction of the lowest 
Landau-level states. 

With Eqs. (|17|l and (|18|l . the operators a± and a|j_ defined in Eqs. H21|l are readily expressed in terms of the original 
coordinates x, y and momenta px,Py. It is not hard to verify explicitly that a± (pooix, y) = 0. The more interesting 
question is the form of the low- lying states ipnoix,y), which are the analogs of the lowest Landau- level states but 
now for an anisotropic rotating trap. In this case, the state ip„Q has n quanta of the plus mode (whose frequency uj^ 
becomes small cx ^/S for (5 = 1 — Q/uJx 0), and zero quanta of the minus mode (which has a finite frequency w_ in 
the same limit). 

The basic relation a+(poo — implies that {Q+/^/2d+)ipoo — —i{d+P+/V2h)Lpoo- Thus a^_^_ipoo = 
{^/2 Q+/d +) if 00 — —i{'\f^d+P^/K)ipQQ, and a straightforward calculation yields 

•^wi^.y) = a\ipQ(i{x,y) = C'/'oo(a;,y), (30) 
where (/Jqo is the normalized ground state, and 



c 



V2{x + iP-y) /2TO7/3+ x + ij3-y 



d+X+{l+ 13+13^) V h l + P+13- 

(31) 
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is a dimensionless complex variable involving a "stretched" combination x + i^-y. As shown in Fig. 1(b), the 
polarization parameter /3_ is real and less that 1 in the limit of rapidly rotating anisotropic trap. This complex 

variable reduces to C = (2; + w)/d'0 for a rotating symmetric trap, where do ~ y^h/{mijJo)- In the more general 
case of a rotating anisotropic trap, the characteristic length that appears in is essentially from Eq. H28|l . 

apart from the common factor + /3_(-/3_) and the additional factor /3_ for y; both these factors remain finite as 

ri — > Wa;. This quasi-isotropic behavior for is very different from the anisotropy seen in the two lengths and Uy 
that determine the x and y structure of the ground-state density |(pooP- 

The higher states within the lowest Landau level have a similar structure. For example, 1^920 = (q^+)^ foo/V^ can 
be written as 

V2o{x,y) = --^ipioix,y) = -^(ifooix^y) 



+ C') </'oo(x,2/). (32) 
The commutator is readily evaluated with Eqs. if^ . ((T7|l . ((T^ and JSIJ, yielding 



which defines the constant c (it depends on the trap frequencies lOx, ujy and the rotation speed SI). In general, 
< c < 1, but it vanishes identically for a symmetric trap since /3± = 1 in this case. Thus (y920 — (C^ ~ c) (^oo/v^, 
namely an even polynomial in C times the complex Gaussian ground state (^oo- As a check on this analysis, note that 
[a+, Q!^)^](poo = V'ooi and direct calculation verifies that [a+,C] = 1- 

The general lowest Landau-level state follows from similar arguments [it is essential here that the commutator 
is a pure number, independent of x and y] 

•^no {x, y) = Pn (C) foo {x , y) , (34) 



where Pn{0 is a polynomial of order n that obeys the symmetry condition Pn(— C) = (~l)"Pii(C)- These polynomials 
are easily obtained recursively with the relation 

Pn+liO^CPniO-C^^, (35) 

with the first few given explicitly as po = 1, pi — C,, p2 = — c, Pa = — 3cC, • • • . The Hermite polynomials Hn(x) 
obey a similar recursion relation j25l | Hn+i{x) — 2x Hn{x) — H'^{x). Direct comparison shows that 

\ TT ( ^ 

r2c 



which readily reproduces the explicit forms given above for small n = 0, • • • ,3. Oktel [I^l] obtained an analogous but 
less general result in the special limit of small anisotropy and rapid rotation. For a symmetric trap (with j3± = 1 and 
c = 0), it follows directly that Pn{C) reduces to the nth power of {x + iy)/do. 

D. Expectation value of single-particle Ho for general lowest Landau- level state 

Let tpLLL — X^n Cn(/3no bc a general linear combination of lowest Landau-level states {ifno}i normalized with the 
condition / dxdy = 1- The expectation value of Hq in Eq. is given by the equivalent operators in l|22|l 

{Ho) = J dxdytpli^i^ ^huj+ (a\_a+ + a+al_^ + ^huj^ (^a''_a- + a^al_^ i^LLL, (37) 

where the angular brackets denote the expectation value with the state ipLLL- Since a-ipLLL vanishes (by construc- 
tion), this quantity reduces to 

(Ho) = \fito^ + \hLo+ J dxdyiplj^j^ (^al_a+ + a+a\_^ tpLLL, (38) 
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where the first term is just the zero-point energy of the unoccupied minus mode. For a symmetric trap, this expectation 
value is readily expressed in terms of the expectation value (x^ + y^) 0, Q- As shown below, a similar but more 
intricate result holds for the rotating anisotropic trap. 

It is convenient to start from Eqs. (|19|) that express x and y in terms of the new canonical variables Q± and P±. In 
turn, these operators are simply linear combinations of the corresponding oscillator variables alj_ and a± , as follows 
from H21(l . For example, 

(a+ + al) f a_ - al) . (39) 



The expectation value of is then given by 

where the cross-terms between plus and minus operators vanish because {a-) = (aL) = 0, and I use the relations 
d\\\. ~ {h/niuj±) [i-o±/ ii±) = h/mjf3±. A similar calculation gives 

and an appropriate linear combination leads to the quantity in Eq. (|38l) . In this way, the desired LLL expectation 
value (Hq) has the simple form 

(Ho) = \hu^ - \hLo+ (^/3+/3_ + ^) + im7c.+ i^p+{x^) + ^{y^)^ . (42) 

For a symmetric trap with lu^ ~ cuy ^ ujq, this result reduces to the familiar LLL expression (Hq) = hCl + 
muo (ujo - fi) (r^) where = -I- j/^ . 

In the case of a symmetric condensate, the special properties of the LLL states yield a simple well-known relation 
between the expectation value of the angular momentum and the mean-squared radius 0, Q 

For an anisotropic condensate, an analogous result follows from the expectation value (L^) = {xpy — ypx) in a LLL 
state. Equations (|19|l for x and y and the corresponding relations for px and py lead to an expression involving the 

combinations + a+a^) and ((at )^ + (a+)^). Comparison with Eqs. (|40() and H41|l and use of Eq. (|15f) then 

yields the following generalization of (|43ll 

For a symmetric trap {lOx = = ^o), it follows by inspection that this result has the correct limit H43|l . because 

7^2tJo, (y') = (^2), and /3_ ^ 1. 



III. INTERACTING GAS IN A ROTATING ANISOTROPIC TRAP 



A dilute interacting Bose-Einstein condensate with N particles in a trap is described by a condensate wave function 
■0 that is here normalized to unity, with J dxdy \ip\'^ = 1. The Gross-Pitaevskii (GP) energy functional for this system 
involves both the noninteracting Hamiltonian Hq from Eq. ^ and the interaction energy 

E[ij] = J dxdy (rHo^J + ^92dN\^\^) , (45) 

where g2d is a two-dimensional coupling constant with dimensions of energy times area. If the condensate is confined 
in a tight axial harmonic trap with oscillator length dz = {rriLOz), then g2d = •s/Stt ?i^as/(mdz), where as is the 
s-wave scattering length 8, 9, 26J. In contrast, for a condensate that is uniform in the z direction with axial length Z, 
the analogous relation is g2d = ATrh^as/ {mZ). The Euler-Lagrange equation for the wave function is the stationary 
GP equation 

i?o0 + 52diV|0|V = MV', (46) 
where the chemical potential /i is fixed by the normalization of tp. 
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A. Lowest Landau-level limit for rapid rotation 



If the trap rotates rapidly, the condensate wave function ijj can be approximated by a general LLL state ipLLL = 
'Y^^CnfnO- In this case, the expectation value of _ffo simplifies considerably to Eq. H42|) . Correspondingly, the total 
energy functional takes the approximate form 



(47) 



2^4^ V 
dxdy 

where I now follow 8] and omit the subscript LLL on the condensate wave function. 

If the restriction to the LLL states is ignored, the absolute minimum of Elll[iP\ is found by varying subject 
solely to the normalization condition. The resulting approximate GP equation 

^^m-yuj+jS+x^ + ^m-yuj+^y^ + g2dN\tlj\'^ = /i (48) 

implies an anisotropic density distribution 

|V,„,„(.,!/)P^^(l-|-^), (49) 

with characteristic condensate radii given by 

Note that the ratio Ry/Rx = P+ oc w+ vanishes as £7 — > Wj,, but the behavior of the individual condensate radii 
requires a study of the chemical potential /u. 



As emphasized by various authors |^ |5|, |a, l3j la ISl j this density is very similar to the familiar Thomas- Fermi form 
for a nonrotating condensate in a stationary two-dimensional trap. In that case, the repulsive interactions expand the 
condensate and reduce the kinetic energy compared to the trap energy and the interaction energy. The situation here is 
very different, because the approximate LLL wave function explicitly incorporates the full single-particle Hamiltonian, 
including the kinetic energy; in this context, the appearance of the squared coordinates x^ and arises from the 
special properties of the LLL states, specifically the result in Eq. H42|l . 

The normalization condition for \'ip\'^ in (|49|l readily yields the condition 



/ m7 u}+g2dN 

^' = \| ^ • (51) 

Note that the chemical potential ^ vanishes proportional to ^/lJ+ oc S^^^ for a rapidly rotating anisotropic trap, where 
(5=1 — n/ujx — > 0. Equation (|49|l then shows that the central density |'0min(0, 0)p also has the same behavior in this 
limit. 

A combination of Eqs. H5()|l and 1)51(1 gives the condensate radii 



Rl - U^^^ - Wj^^. (52) 



Correspondingly, the normalized minimizing density H49|l becomes 



IV'min(a:, y)p = -2— f 1 _ ^ _ ^ ) . (53) 



nRxRy \ Rx Ry 



Since uj+ and /3+ are both proportional to \/S as (5 ^ for fixed trap anisotropy, it is clear that R^ grows like S '^1^ 
and i?y shrinks like 5^1^ for small 5, which reflects the conservation of total number of particles. In particular, the 

total area irRxRy diverges like i-^^^^^ oc (J^^/**. This anisotropy of the minimizing TV-body condensate in the limit 
(5 — > is quite different from the anisotropy of the LLL single-particle ground state ipoo, where Eq. H28|l shows that Ux 
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grows but ay approaches a constant as 5 — > 0. For the minimizing density (|53|l . it is straightforward to evaluate the 
mean-squared displacements {x^) — and (y^ 



\B?y The mean angular momentum in Eq. 1441) then becomes 



/3+ 



J- 



Ri 



/3- 



(54) 



The condition .g2d"-(0) <C for the validity of the lowest Landau-level approximation can now be made explicit. 
Use of 52d?T-(0) ~ 2g2dN I {iiRxRy) from Eq. for the condensate radii and the relation g2d ~ ^TrTi^as/ {mZ) 
for a uniform condensate of thickness Z yields 



Nas 
Z 



< 1. 



(55) 



Since uj^ — > and w_ remains nonzero for sufficiently rapid rotation ^ ujx), this condition can always be satisfied. 
In the special case of a symmetric trap, the minimizing density has the isotropic form Q 



|V'min(r)|' 



r 



(56) 



with 



i?o - 



-,1/2 



(57) 



This squared condensate radius diverts as 17 - 
(|47|l for a symmetric trap becomes Q ^^lllI 



limit. 



/dQ — h (ii?§/c?Q — l); this result agrees with Eq. because 



LOo. Similarly, the absolute minimum of the LLL energy functional 
= hfl + ^itlujqRI. As a simple check, it is easy to verify that 



^Rq in the symmetric 



B. Density of vortices 

The general LLL state ipLLL is a linear combination of the states (pno- Apart from a normalization factor l/VrJ., 
each of these states is a polynomial Pn(C) multiphed by the ground state ipoo, where ( from Eq. is proportional 
to X -|- iP-y. Thus the general LLL state also involves a polynomial in ^ that can formally be factorized to write 

'ipLLL oc (fooYliC - Q) ■ (58) 
i 

The corresponding LLL particle density riLLL = IV'lllP becomes 

Apart from an additive constant, the logarithm of this relation gives P.l8ll2^l2^ 

,2 



y^lii\{x - Xjf + (3^_{y ~ yj'f\ = + IL + lnnLLL{r) 



1^2, P~y 



i2 



+ 



x^ + '—^\+\nnLLL{r). (60) 



Here I use Eqs. (|28|) and (|29|l for the anisotropic ground state and note that — P+. 

To include the anisotropy of the complex variable (x x + i(3^y, it is convenient to introduce the rescaled variables 
x' = X and y' = f3-y- Application of the rescaled Laplacian V'^ = d^/dx'^ + /dy'^ readily gives 

^ V- In |r' - r-^f = ^ (l + ^) + 1- -"^(-)- (61) 
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Since V'^ln|r'p = An S{x')S{y') = {An/ f3^) S{x)S{y), this relation implies 

nAr)^^ + ^VHnnLLL{r), (62) 

where ny{r) = S^'^^{r — rj) is the two-dimensional vortex density, and I have again used Eq. H28|) . For a rapid 
rotation speed ^ ^ u>x ^ ujy and any trap anisotropy, the frequency 7 = (w^ — w^)/(2f2) is given in terms of 
elementary expressions from Eqs. l|2Jl and Q), as shown in Fig. 1(a) for a typical anisotropy ratio LOy/uJx — 1-2. 

To estimate the vortex density in Eq. (|62(l . assume that the equilibrium particle density ulll is that given by 
the absolute minimum solution in Eq. (|49l) . with nLLL{r) (x 1 — x^/R^ — y^/Ry. A straightforward calculation then 
yields 



ni'-f (3. 
nv[x,y) « — ^ 



27rh 2n{l-xyRl-yyRl) 



1 1 



Rl PlRl 

(63) 



For an isotropic trap {lOx — = ^0)7 this expression reduces to the well-known axisymmetric result [7l l27il2^ 

"^^^^"^-i!(i-.W' ^''^ 

where i?Q is given in Eq. H57I) . 

Comparison of these two expressions shows some interesting differences: 

1. For fixed angular velocity < wo, the vortex density 1)64(1 in a rapidly rotating symmetric trap decreases 
gradually and isotropically away from the center of the condensate. Such behavior has been observed at lower 
angular velocities in the mean-field Thomas- Fermi regime |29l |. In contrast, the general expression for the vortex 
density in Eq. (|63|l displays explicit anisotropy between x and y. 

2. To sharpen this analysis, it is convenient to focus on the central vortex density 



21:% 27r V^x P-^l 

For a symmetric trap, Eq. H64() shows that n„(0) increases monotonically with increasing f2 < cjo, because the 
mean condensate radius _Ro grows in the same limit. For an anisotropic trap, in contrast, the frequency 7 and 
R~^ both decrease with increasing Q,, whereas increases. Thus the combined effect of anisotropy and rapid 
rotation can, in principle, yield a central vortex density nt,(0) that varies non- monotonically with increasing O, 
as can be seen for typical numerical examples. 

It is not clear whether either of these behaviors would be observable in practice. 

The present discussion has focused on the "macroscopic" parabolic density profile that provides an absolute mini- 
mum of the energy in the rotating frame, ignoring the local distortions associated with the vortex cores. In practice, 
these phenomena have very different length scales: in the rapidly rotating limit, the vortex core and the intervortex 
separation are both of order ^J'h/ {mfl), whereas the condensate radii Rx and Ry are generally much larger. Thus it is 
possible to treat the condensate density as locally uniform over the size of an individual vortex. If the vortex lattice 
is treated as triangular and unbounded, this analysis yields a simple renormalization 0, S IS of the interaction 
constant g2d — > bg2d, where b « 1.1596 is the numerical value for a triangular Abrikosov vortex lattice (SliJ . Apart 
from this rescaling of the interaction parameter, the description remains essentially unchanged. 



IV. CONCLUSIONS AND DISCUSSION 



This work has examined the behavior of a two-dimensional Bose-Einstein condensate in an anisotropic harmonic 
trap (with general trap frequencies tOx < Wj.) that rotates rapidly at an angular velocity Q ^ lOx- The single-particle 
Hamiltonian in Eq. is exactly soluble [wL iTj . Il3| , although the detailed form of the low- lying quantum- mechanical 
states in Eq. (30) for an arbitrary anisotropy has apparently not appeared previously. 
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The ground-state wave function ipQQ{x,y) in Eq. H27() is an anisotropic Gaussian with fi-dependent characteristic 
lengths and Uy given in (j^SJ. In addition, the ground state has a phase proportional to xy [2^ In the 

quantum problem, this behavior reflects the classical velocity potential for an irrotational (vortex-free) fluid confined 
in a rotating elliptical boundary . 

Similar to the case of a rapidly rotating symmetric condensate, the two eigenfrequencies u!± in Eq. (0) for a rapidly 
rotating anisotropic condensate have very different magnitudes: LU-f- vanishes as lUx (the smaller of the two trap 
frequencies), but lu^ remains finite in this same limit. If the mean interaction energy ~ g2dn{0) is small compared to 
the gap hu- between the ground state and the first excited Landau level, then the system can be described with the 
set of lowest Landau-level states ipno{x,y), where ri is a non-negative integer describing the number of plus quanta, 
each with the small energy Tilo^. Apart from the Gaussian factor (^00(2^1 ?/)j these states involve a polynomial p„(C) 
of order n, where C, single complex variable proportional to a; + iP-y and /3_ < 1 is a positive constant for any 

Within the set of lowest Landau-level states ipLLL — J^n'^ri'PnQ, the expectation value (Hq) of the single-particle 
Hamiltonian can be reduced to an anisotropic linear combination of (a;^) and (y^). An anisotropic Thomas-Fermi-like 
density |V'min(a;, ?/)P oc 1 — /R^ — y^ / provides an absolute lower bound for the total energy of the interacting 
system. Here the condensate radii Rx and Ry are given in Eq. ()52|l : R^ grows and i?^ shrinks as — > w^;, while the 
area irRxRy of the elliptical condensate grows slowly in the same limit. 

As emphasized by Ho the particle density nLLL{x,y) = \ipLLL{x,y)\'^ in the lowest Landau-level limit also 
contains a description of the associated vortex density nv{x,y). This situation arises because the general linear 
combination of lowest Landau-level states is essentially a polynomial P{C) iu the single complex variable C oc x + i(3-y, 
and the zeros {Q} of P{C) represent the positions of the vortices in the xy plane. The actual particle density ulll 
has small-scale structure arising from the vortex cores, superposed on the parabolic global shape. If this fine-grain 
aspect is ignored, the resulting vortex density follows immediately in Eq. Ht)3|l . 

This work raises several interesting questions: 

1. Feynman's familiar expression for the mean vortex density np = mfl/^nh) in a large symmetric rotating 
condensate requires modification for an anisotropic rotating condensate because of the irrotational fiow induced 
by the rotating walls. Specifically, this irrotational flow contributes to the total angular momentum and thus 
lowers the energy E' = E — JlLj in the rotating frame, delaying the flrst transition to a state containing a 
single quantized vortex |2^. Indeed, the measured critical angular velocity flc for the appearance of the first 
vortex in a rigid elliptical or rectangular cylinder containing uniform superfiuid ^He exceeds that for a circular 
cylinder by a factor that increases with increasing anisotropy ,32i| . confirming the theoretical prediction. A 
similar but more complicated situation occurs for slowly rotating weakly interacting anisotropic Bose-Einstein 
condensates [llj . Although the analogous situation with many vortices has not been studied in detail, the 
increased critical angular velocity for the appearance of the first vortex suggests that the mean vortex density 
at moderate rotation speeds in an anisotropic trap is likely to be smaller than in a corresponding symmetric 
trap. 

In contrast, for a rapidly rotating anisotropic condensate, Eq. (|65|l gives the central density ny{0) « m'y / {2ttTi) , 
apart from finite-size corrections associated with the condensate radii. This value of ny{Q) typically exceeds np: 
because 7/(2ri) — > {ioJ^ + Wy)/(4(jj^) > 1 in the limit Vl ^ ujx- The physical basis for such enhanced vortex 
density is not immediately obvious, and additional clarification would be desirable. 

2. The present LLL approach ignores the fine-grain structure of the vortex lattice (apart from the renormalization 
of the interaction parameter g2d by the "Abrikosov parameter" b ~ 1 . 1 6) 0, 0, H, 0, 113 • In the special case of a 
nearly symmetric trap rotating close to the limit of instability, Oktel 12] finds that the vortex lattice remains 
almost exactly triangular, with the principal lattice planes aligned with the direction of weak trap confinement. 
It is not clear whether this situation holds for arbitrary anisotropy and less extreme rotation speeds. 

3. As f2 approaches the weak confining frequency lOx, the condensate becomes essentially one dimensional and the 
vortices must then rearrange themselves to form a one-dimensional structure "23,123]. Ultimately, this behavior 
may be pre-empted by some sort of transition to a correlated (nonsuperfluid) state, as has been predicted for a 
symmetric trap |33| . 

It will be interesting to investigate these various questions in detail. 
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